Geometric phases of d-wave vortices in a model of lattice fermions 
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We study the local and topological features of Berry phases associated with the adiabatic transport 
of vortices in a d-wave superconductor of lattice fermions. At half filling, where the local Berry 
curvature must vanish due to symmetries, the phase associated with the exchange of two vortices is 
found to vanish as well, implying that vortices behave as bosons. Away from half filling, and in the 
limit where the magnetic length is large compared to the lattice constant, the local Berry curvature 
gives rise to an intricate flux pattern within the large magnetic unit cell. This renders the Berry 
phase associated with an exchange of two vortices highly path dependent. However, it is shown that 
"statistical" fluxes attached to the vortex positions are still absent. Despite the complicated profile 
of the Berry curvature away from half filling, we show that the average flux density associated with 
this curvature is tied to the average particle density. This is familiar from dual theories of bosonic 
systems, even though in the present case, the underlying particles are fermions. 



PACS numbers: PACS 

Introduction. A phenomenology based on a BCS-like 
pairing state with d-wave symmetry has led to consid- 
erable success in understanding the properties of quasi- 
particles in high-T c superconducting cuprates. This in- 
cludes the mixed state of these systems, where a magnetic 
field H Cl < H < H C2 is applied and leads to the presence 
of an Abrikosov vortex lattice. Effective models^"^ have 
been developed that describe the dynamics of the quasi- 
particles under the simultaneous influence of magnetic 
field, the supercurrent flow due to the vortices, and in 
some cases the underlying microscopic lattice. The vor- 
tices of the mixed state are usually assumed to be static, 
i.e., frozen into the Abrikosov lattice. However, it has 
been argued^ that both as a result of the small coher- 
ence length £, and possibly the proximity to an insulat- 
ing state, fluctuations of vortices may play a fundamental 
role. Moreover, at sufficiently high magnetic fields below 
H C2 , it has been predicted that thermal and/or quantum 
fluctuations may melt the vortex lattice or glass, leading 
to a "vortex liquid" regime^. For all these reasons, it is 
desirable to construct effective theories that include the 
vortices as fundamental dynamical degrees of freedomPD 
Such a construction is readily available in systems where 
the constituent particles are bosons, through the well- 
known Kramers- Wannier duality 19 . In Fermi systems, 
however, the vortex degrees of freedom only exist as 
dual partners of bosonic Cooper pairs that are them- 
selves emergent particles. This arguably complicates the 
task of passing directly from a microscopic description 
in terms of electrons to an effective theory in terms of 
vortices, requiring more ad hoc assumptions. Such effec- 
tive theories have been previously discussed in a contin- 
uum formalism.^ In this paper, we aim to establish some 
key parameters of these theories in a microscopic lattice 
model. This is similar in spirit, but physically different, 
from earlier considerations for bosons in the absence of 



a lattice--^. These defining parameters include the quan- 
tum curvature felt by the vortices in the condensate, that 
is, the effective magnetic field experienced by them, and 
their mutual statistics. Specifically in the d-wave pair- 
ing case, where the continuum description of vortices is 
somewhat plagued by subtleties concerning self-adjoint 
extensions^, our microscopic starting point also serves 
as lattice regularizatiorP, which allows for a controlled 
study of the desired universal properties. 

Model description. We will study the Berry phases^ 
of vortices in the BCS-Hofstadter model, which has been 
used previously as a microscopic description of the mixed 
state in the cuprates^. 
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In the above, the sum (rr') is over the nearest neigh- 
bors, and the hopping terms are just those of the Hof- 
stadter model, to be described below. The correspond- 
ing uniform magnetic flux through the plaquettes of the 
lattice mirrors the fact that the penetration depth is 
much larger than the coherence length, as befits a type 
II superconductor. Assuming symmetric gauge A(r) = 
(— y/2, x/2)&, where $ the magnetic flux per plaquette, 
the hopping amplitude assumes the form 

t rr , =te- i K' Adl , (2) 

where r refers to the discrete sites of the lattice. The 
d-wave paring term is defined as 

Arr' =r/ r _ r <A arr ,e je "' (3) 
W = +/- if (r - O || x/y (4) 

Here, ^ r _ r ' encodes the d-wave symmetry. Ao, rr ' is es- 
sentially constant, except for a suppression of amplitude 
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near the vortex core. We follow RefP in defining the 
pairing phase factor e 6rr ' via 



| e i<£(r) 4- e i4>(r>) | ' 
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i.e., as a link-centered average of a field <j>{r) that satisfies 
the following continuum equations, 

V x V0(r) = 27r£^(5(r- r t ) (6a) 

i 

V • V</>(r) = , (6b) 

where the denote the vortex positions, which can take 
on continuous values. The total number of vortices ny 
is equal to the number of half flux-quanta $o> lxly& — 
ny$ where <& = 7r in natural units, and l x , l y are the 
number of unit cells in the x and y direction, respectively. 
Eq. ffj} can be solvecP via 

<K r ) = ^2{arg[(r(z-Zi, u, u / )]+2-y(x-x l )(y-y t )+v Q -r} , 

i 
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where cr(z, ui, u>') is the Weierstrass sigma-function with 
half periods u) — l x /2, uj' — il y /2, z = x + iy, and 
the sum is over vortex positions. Integration constants 
v = 2 J2i A(r -) and 7 = ^ — f- have been chosen such 
that the superfluid velocity V5 = V(j>/2 — A satisfies pe- 
riodic boundary conditions, and averages to zero over the 
magnetic unit celP, and 77 = C( w ) is pure imaginary, with 
C the Weierstrass zeta-function. 

The pairing phase factor e l6rr ' in Eq.fJHJ is ill-defined 
when the denominator goes to 0. This is unacceptable 
since we mean to continuously change vortex coordinates 
in the following. To remove this singularity, we define 
Ao,rr' as 



A , rr ' = A [l - exp(-- 



Calculation of the Berry phase. In the following, we 
will consider the model Eqs. ([T])-([9) as a function of vor- 
tex positions {r,;}. We note that the simultaneous pres- 
ence of the magnetic field and the discrete ionic lattice 
generically opens up a gap in the quasi-particle spectrum 
of the (i-wave superconductor, except for special vortex 
configurations that respect inversion symmetr . The 
Berry phase associated with the motion of vortices is 
thus well-defined. We further remark that the model 
defined above is traditionally studied by means of a sin- 
gular gauge transformation 2 , that, on average, removes 
the magnetic field. This is inconvenient for present pur- 
poses, since the precise transformation depends on vor- 
tex positions, and the Berry phase is clearly not invariant 
under unitary transformations that vary along the par- 
ticular path in question. We thus need to stay within 
the present framework of magnetic translations and as- 
sociated boundary conditions. 

To study the Berry phases associated with the motion 
of vortices, we first note that within our model the vor- 
tex positions are well-defined continuous parameters that 
are, at least for large enough lattice, entirely encoded in 
the pairing amplitudes A rr < . The Berry phase associated 
with vortex motion along closed paths may be computed 
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where the |f2j) are the ground states of the system along 
a reasonably fine discretization of the path. The above 
formula has the advantage (over the standard integral 
formula) that a random, discontinuous phase that each 
acquires in numerical diagonalization automatically 
cancels. Each ground state is constructed as the vacuum 
of Bogoliubov operators 
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where Ao and £ are constant parameters. This leads to a 
suppression of pairing amplitude on links near the vortex, 
hence £ may be thought of as a core radius. 

We further impose periodic magnetic boundary condi- 
tions on our model as follows: 



Re 
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(9) 
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In the above, the magnetic translation operators T x and 
T y are defined by letting R = xory. We note that with 
the boundary conditions ([9| imposed on electron opera- 
tors, the physics is also periodic in the vortex positions 
rj. That is, one may see that the formal replacements 
r-j — » Yi + l x x, ri — > Yi + l y y affect the Hamiltonian by 
a unitary transformation, as given explicitly below. In 
particular, the quasi-particle spectrum is invariant under 
such replacements. 



where the matrices U rn — ii n (r), V rn = v n (r), satisfy 
Bogoliubov-deGennes equations 

-A-" ,t„) i-v) - i-v) <"> 



for non-negative eigenvalues E n . It is clear from Eq. ( 12 ) 



that the state |0) = n r c ril0) is a vacuum of both the 



operators 7„^ and 7^ where |0) is the vacuum of the 
c rcr operators. The ground state of the Hamiltonian thus 
can be constructed as 



N>=n^4io) 



(15) 



Using this last relation, and the inverse of Eq. (12), one 
readily obtains 



(n i \n j ) = det{u i uj + ViVj). 



(16) 
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Results. We first consider the important special case 
of Eq. ([!]) with fi = 0, or half-filling. In this case the 
Hamiltonian is invariant under the anti-unitary charge 
conjugation operator defined via Cc ra C = (— l) r cJ CT , and 
the unique ground state is then invariant under C 
as well (up to a phase that c an be made trivial). It 
then follows directly from Eq. (11) that e n = ±1. The 



first immediate conclusion from this is that as long as 
vortices are moved along contractible paths, the Berry 
phase must be +1 for continuity reasons. If vortices were 
hard-core particles this would, in principle, still leave the 
possibility of fermionic statistics. However, careful exam- 
ination shows that the Hamiltonian can be analytically 
continued without difficulty into configurations were two 
vortices fuse into a double vortex at a given location. 
Exchange paths are thus contractible, and hence vortices 
must satisfy bosonic statistics. We have tested this for 
various lattice sizes and exchange paths. The model does, 
however, become singular when vortex positions are for- 
mally approaching lattice sites, see Eq. Q. It is thus 
possible that lattice sites carry an effective 7r-flux felt by 
vortices encircling such sites. We have carefully checked 
that this is not the case in our model. Hence at half fill- 
ing, all Berry phases are unity. The above observations 
also hold for the s-wave case. 

The observation that vortices are bosonic is non-trivial, 
since time reversal symmetry is absent, and hence gener- 
ically in two spatial dimensions even non-Abelian statis- 
tics are possible, as is the case if the pairing symmetry 
is p + ip22i. Indeed, when we move away from half filling, 
there is no longer any symmetry that requires the Berry 
phase to be trivial. We will now show that this situation 
leads to a very intricate landscape of non-trivial quantum 
curvature. 

The Berry curvature is defined as the Berry phase 
around an infinitesimal area, divided by the size of this 
area. In the following, we consider a lattice containing 
only two vortices in the presence of periodic boundary 
condition. One vortex remains fixed, while for any point 
within the unit cell, we calculate the Berry curvature as- 
sociated with the motion of the other vortex according 
to Eq. ( 11 ). The Berry phase around arbitrary loops can 



be obtained as the integral of the Berry curvature over 
the enclosed area. The result for a 12 x 10 lattice at 
/i = 0.05 is presented in Fig. [lja). It is apparent that 
the Berry curvature in this model is a highly non-trivial 
function of position for any /j, ^ 0. One observes that the 
curvature is conspicuously concentrated on the links and 
the sites of the lattice, even though the vortex positions 
themselves are formally not tied to the discrete lattice. 
Singular structures form in particular around the lattice 
sites. These are described by B(r) ~ ai5(r—Ti)+fi(6)/r, 
where B(r) is the curvature, and 9, r refer to polar co- 
ordinates with the lattice site rj at the origin. The pa- 
rameters a, and the functions fi depend sensitively on 
details such as the lattice size, /x, the site index i, and 
the position of the other, fixed vortex. Yet another in- 
teresting feature is the structure seen in the vicinity of 
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FIG. 1. (Color online) Berry curvature for 12 by 10 lattice 
in the presence of two vortices, for [i — 0.05. The calculated 
curvature at a given location is associated with the motion 
of one of the vortices at this location, while the position of 
other vortex is held fixed, (a): 3D view of the Berry curva- 
ture in the vicinity of one plaquette. Lattice sites and links 
are clearly visible, since the curvature is chiefly concentrated 
on the latter, with singularities at the sites. The remaining 
space within the plaquettes is nearly flat, (b) Top view of 
the lattice. The position of the fixed vortex is noted by an 
asterisk. No singularity is present at this position. However, 
an interesting structure reminiscent of the shape of a d x 2_ y 2 
orbital appears in the vicinity, (c): 3D view of the entire 
lattice. Detailed features appear highly irregular, except for 
the fact that the curvature is everywhere concentrated on the 
sites and the links of the lattice, and on the structure seen 
around the position of the fixed vortex. 



the fixed vortex, which is somewhat reminiscent of the 
shape of a d x i_ y i orbital. However, this structure does 
not seem to be reflective of by the pairing symmetry, 
but rather more the lattice symmetry, as similar calcu- 
lations for the s-wave case show. We note that again no 
singularity indicative of a flux tube carried by the fixed 
vortex appears in Fig.JlJb) at the position of this vortex. 
This implies that we should still think of these vortices as 
bosons, which move in an effective background magnetic 
field. 

The complex nature of these features and the strong 
sensitivity on model parameters are likely yet another 
facet of the fractal nature of the physics of the Hofstadter 
model. To wit, in view of the fractal nature of the wave- 
vector dependence of spectral features of the Hofstadter 
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model, it is reasonable to expect that the response to a 
spatially inhomogeneous perturbation (coupling to many 
different wave vectors) is characterized by complicated 
and possibly chaotic spatial modulations. The addition 
of a pairing order parameter with vortices clearly repre- 
sents such a perturbation. Here we are mostly interested 
in how to reconcile the complex features seen at fi =/= 
with the trivial ones seen at fi — 0. It is clear that our 
ability to precisely define the vortex position on scales 
below the lattice constant is dependent on conventions, 
even though in the present case a natural convention is 
available, since our ground states are naturally param- 
eterized by the vortex positions in the continuum field 
^ used to define the Hamiltonian. We have, however, 
tested the robustness of the qualitative features shown 
in Fig. [I] by varying the precise form of the pairing or- 
der parameter Eq. In particular, we have varied the 
core parameter £ in Eq. ([8]), and tested various alterna- 
tive forms for Eq. We also introduced variations in 
the boundary conditions described above. In all cases we 
found that the qualitative features of the Berry curva- 
ture remained unaltered. Although we believe that the 
curvature landscape of Fig. |ljc) is interesting in its own 
right, it is appropriate to make this landscape subject 
to some coarse graining procedure. It is interesting to 
ask whether such coarse graining leads to a recovery of 
one of the basic facts suggested by conventional wisdom 
about vortex-boson duality, namely that the curvature 
discussed above is directly tied to particle (here, Cooper- 
pair) density. We will show in the following that this 
statement is recovered when the Berry curvature is av- 
eraged over the magnetic unit cell (as opposed to the, 
typically much smaller, lattice unit cell). To this end, 
we again consider a lattice containing only two vortices 
within a single magnetic unit cell, subject to the bound- 
ary conditions ([9|. 

Let the coordinates of the "moving" vortex be r = 
(x, y). As remarked initially, a formal shift of x by l x 
changes the Hamiltonian by a gauge transformation. We 
have H -> W X HU X with U x = e ^(i-y/i v )N/2 ^ where ft ig 
the particle number operator. Analogous relations hold 
for y -> y + l y , with U y = e Ml+*/**)#/2, We nQW calcu _ 
late the Berry phase associated with a rectangular path 
of dimensions l x , l y around the lattice. We may then 
choose a ground state phase convention |0(r)) along the 
path satisfying 



\Q(r + l y y)) = Ul\Q(r)) 

\n(r + i x x)) = ul\n(v)) 



(17) 



along the horizontal and vertical path segments, respec- 
tively. The consistency of Eq. ( 17 ) with the continuity of 



|£l(r)) along the path follows from the observation that 
for any ground state, UlU*U y U x \n) = The lat- 

ter holds because U^U^.U y U x — exp(inN) and because 
the ground state Eq. ( 15 ) always has even particle num- 
ber parity. We determine the Berry phase for the rect- 
angular path as the integral over the Berry connection, 



(f2(r)|V|f2(r)), where we observe that 

(fi(r + ka)|V|n(r + l x x)) = (fl(r)|l/ x VE£|fi(r)) 

= <fi(r)|V|n(r)> + ^<JV)0 



(18) 



It is clear that only the last term survives a cancella- 
tion between the vertical path segments, giving iw(N) /2. 
The same contribution is obtained from the horizontal 
segment. We thus obtain 7 = ir(N), or 2ir times the 
number of Cooper pairs in the system, in agreement with 
general expectations based on duality arguments applied 
to Cooper paira^. It is worth noting that the quantity 
7, when expressed as an integral of the Berry curvature 
over the entire lattice, is formally reminiscent of a Chern 
number. It is not truly a Chern number, though, since 
the boundary conditions (17) do not quite allow one to 



make contact with one-dimensional vector bundles over 
the torus. Indeed, 7 is not quantized, as (N) may take on 
arbitrary values in [0, 2l x l y ]. We note that the derivation 
above is independent of the pairing symmetry. 

Conclusion. The present study establishes several as- 
pects of Berry phases associated with vortex motion in 
a microscopic model of superconducting lattice fermions. 
It is shown that these vortices behave as bosons which, 
away from half filling, are subject to a non-trivial ef- 
fective magnetic field. In an average sense, it has been 
shown that this effective field is tied to the density of 
Cooper pairs. This is expected based on boson-vortex 
duality, and was seen to emerge here in a microscopic 
model of fermions. We emphasize that the simple re- 
lation between Cooper-pair density and effective field is 
only seen to emerge after averaging over a magnetic unit 
cell. This may be used to justify a direct proportionality 
between Cooper pair density and Berry curvature in the 
long wavelength effective theory. However, our results 
also indicate that care must be used in order to justify 
such a relationship in general. On the one hand, this is 
true because of the relatively large non-uniformity of the 
observed Berry curvature within the magnetic unit cell. 
Moreover, in the presence of particle hole symmetry we 
have found that the Berry phase associated with closed 
paths is always zero, and thus corresponds to it times 
the average enclosed particle number only for such paths 
that happen to enclose an even number of lattice sites. 
In this case, the background field appearing in the effec- 
tive theory should clearly be zero, and should not follow 
the total Cooper pair density. This result will be robust 
to small perturbations respecting particle hole symmetry, 
and is thus true for a wide class of microscopic models. 
We conjecture that the complex landscape of the Berry 
curvature away from half filling is a facet of the fractal 
properties of the Hofstadtcr model, and believe that it is 
worthy of further investigation. 
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